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Abstract 



> 

^ | ' We study the tunneling through an oscillating delta barrier. Using time peri- 



odicity of the model, the time-dependent Schrodinger equation is reduced to 



£NJ ' a simple but infinite matrix equation. Employing Toeplitz matrices methods, 

' the infinite matrix is replaces by a 3x 3 matrix, allowing an analytical solution. 

o 

' Looking at the frequency dependence of the transmission amplitudes, one ob- 

serves a new time scale which dominates the tunneling dynamics. This time 
scale differs from the one previously introduced by Biittiker and Landauer. 
> ■ The relation between these two is discussed. 

c 

us 
cr. 

> : 

• i-H , 

^Xj ■ I. INTRODUCTION 

.5?; 

The question regarding the actual time spent during tunneling through a barrier has 
been in discussion for over half a century |L]-§|. The introduction of high-speed tunneling- 
based semiconductor devices in recent years has brought new urgency to the problem. 
However, in spite of the long history of this subject it still remains controversial. 

Most of the opinions on this matter can be classified into three categories. The most 
apparent way to define and calculate tunneling times is through the dynamics of the wave- 
packet. Some suggestions use the stationary phase approximation and obtained the time 
the peak of the wave packet spent in the barrier region 0]. In a more recent versions of this 
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approach, the peak is replaced by the centroid of the packet ||. Nevertheless, this approach 
is problematic. It has been shown || by an explicit example that a tunneling situation can 
be set up such that the peak or centroid of the transmitted packed emerges before the peak 
or centroid of the incident packet has even arrived to the barrier, leading to negative time 
and demonstrating the lack of causal relationship. 

A more sophisticated way involves the introduction of a physical clock which is used to 
determine the time elapsed during tunneling. One can look at either the effect of the clock 
on the tunneling particle or at the clock variable itself. Mainly, three versions of clocks 
have been proposed and investigated in the literature. The first clock approach studies the 
precession of the spin of the tunneling particle due to a uniform infinitesimal magnetic field 



confined to the barrier region [10|. It was later noted that in addition to the precession in 



the plane perpendicular to the magnetic field, the spin-up component in the magnetic field 



direction tunnels preferentially due to the Zeeman splitting [|TTJ. These two processes caused 
by the field can be used as physical clocks. One gets two different time scales from these 
two. The two time scales are equal to the real and imaginary parts of the complex time 
obtained using the Feynman approach [T^]. A second clock approach uses an incident wave 
composed of two interfering waves [|13|| . If these two are of opposite spin direction, the total 
incident particle flux is uniform in time with an oscillatory spin. 

Biittiker and Landauer (BL) studied tunneling through a rectangular barrier with a 
small oscillating component added to the height ||. The incident particles with energy 
E can absorb or emit modulation quanta hQ during tunneling, leading to the appearance 
of sidebands with energies E ± nMl and corresponding intensities I± n . Looking at the 
Q dependence of these quantities, BL determined the critical frequency as that in which 
the transition between the behavior at the two limits of quasi-static barrier and the average 
barrier occurs. The inverse of this frequency is an indication of the time-scale of the tunneling 
process. 



Third, there are many who object the question in the first place It is often said that 



the concept of tunneling time is not well-defined in the context of quantum mechanics, and 



does not correspond to any observable. This is related to the absence of a time-operator 
in quantum mechanics. This problem can be bypassed using an operator which measures 
whether the particle is in the barrier or not, and then (in time-dependent problems) aver- 
aging over time. The result, divided by the incident flux, is termed the "dwell time" |TT 



Nevertheless, this solution is far from being satisfactory, since taking expectation values, 
one can not distinguish between the transmitted and reflected particles, and thus the time 
obtained is just an average over two (possibly different) times related to these two processes. 
This problem follows from the absence of a well-defined history of a quantum particle |15 



As far as the particle was not measured in the barrier, it was not there at all and thus 



formally the time spent in the barrier is zero [16 



In this paper we do not want to take a stand in this debate, but rather to point out 
that even if one argues that the actual time spent inside the barrier region is not defined or 
zero, it is agreed that the wave function spent a finite time inside the barrier. The question 
about the time-scales related to the wave-function transition through the barrier is at least as 
important as the (maybe more fundamental) former one. After all, the quantum dynamics is 
determined by the wave-function and thus many questions in which time scales are relevant, 
are really not related to the actual time as measured experimentally, but rather at the time 
in which the wave-packet - describing the particle in times it is not measured - was in the 
regime in interest. This approach is motivated by the physical clocks approach. Again, even 
if the application of the proposed experiments to the 'strong' question of tunneling time is 
arguable, one can't deny the need for tunneling time scales while considering the process 
in a time-dependent environment such as considered in the BL approach. In what follows 
we use the same approach and consider the tunneling process through an oscillating delta 
barrier 

V(x,t) = V 6(x)(l + ecoa(Qt)). (1) 

As in the BL case, the tunneling particle can absorb or emit any number of energy quanta 
HQ while being transmitted or reflected (see Fig. |l|). Using the time periodicity of the 



Hamiltonian, the Schrodinger equation is reduced to the inversion of an infinite tridiagonal 
matrix. In practice, the matrix can be truncated very efficiently and reliable numerical 
solutions are easily obtained. In the large quantum numbers regime the three-diagonals 



elements become almost constant, and the matrix is approximately a Toeplitz matrix |]IT,IB 



In this case one can reduce the problem to the inversion of a 3 x 3 matrix, and obtain 
analytical expressions. 

Following Ref . [p|,[i3| , we study the frequency dependence of the transmission amplitudes 
in the limit of small oscillation amplitude e. We show that there exists a time scale r, such 
that the side-band asymmetry (studied by BL) is a function of Qr. This time scale is not 
the same as the one defined by BL which vanishes in the delta barrier limit. This indicates 
the possibility of more than one time-scale for the tunneling process. 



II. MODEL 
A. Analytical Solution 

We consider a time periodic potential (fl|) where Vq is the barrier strength, e (0 < e < 1) 
is the modulation strength amplitude and Q the barrier frequency. 

We first note that since the problem is periodic the general solution of the wave equation 
is a superposition of eigenstates of the Floquet operator, i.e., the operator which shifts the 
solution in time by one period. Thus we can solve for each such eigenstate separately. Since 
the frequency of the incoming part of the asymptotic behavior of the solution is determined 
by the incident energy, this energy fixes the corresponding eigenvalue to be e - lEiT l h , where 
Ei stands for the incident energy. The outgoing parts of the asymptotic solution correspond 
to the transmitted and reflected particles, and their frequencies must match the same eigen- 
value, i.e., £- %E f T l n = e -iEiT/h^ wnere £j j s f| na l energy, after transmission or reflection. 
We thus conclude that the frequency of all the asymptotic outgoing waves should differ from 
the initial frequency by an integer multiple of Q, i.e., the energy change is an integer number 
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of energy quanta HQ. 

The corresponding Schrodinger equation is 

On both sides of the barrier the solution of Eq. (0) is just a free particle solution. We thus 
separate the solution, and write it on both sides as a superposition of free particle functions, 
satisfying the scattering boundary conditions, and having appropriate frequencies: 

ipi(x,t) = e {ik ° x - iuJot) + J2n r ne { ~ iknX ~ iuJnt) x < 
iP(x,t) = { ' (3) 

1p r (x, t) = En t n e^ knX ~ luJ ^ X > 

where r n and t n are the reflection and transmission coefficients, uj n = (E + nhQ)/h, 



k n = \j2muj n /% and the summations are over every integer n. For negative uj n , k n become 
imaginary, k n = iK n . The corresponding wave-functions can not be extended to the whole 
space, but still can contribute to the wave-function on one side. One picks the decaying 
solution to ensure integrability. 

Continuity of the wave function requires i/ji(0,t) = i[) r (0,t), implying the following rela- 
tions, 

r n = t n , n^O 

(4) 

1 + r = to, n = 0. 

Integrating Eq. (Q) over the barrier region one obtains a tridiagonal matrix equation : 

£ £ 

(2ik n /B - l)r n - -(r n+1 + r n _i) = <5 , n + ^(So, n +i + S , n -i), (5) 
where B = 2mVo/h 2 . 

In order to solve Eq. (|^) we need to invert an infinite matrix. Practically, it is enough 
to truncate the matrix after about a dozen rows around the center (n = 0), since the far 
r n 's approach zero quickly (we denote this solution by FS, for full solution). Moreover, it 
is easy to show that for |n| 3> 1 the matrix elements are approximately constant. In this 
case we can treat our matrix as a Toeplitz matrix |l7j , and assume that the solution decays 



exponentially. The solution is then obtained by inversion of a 3 x 3 matrix (we call this 
solution TS for Toeplitz solution). 

Let us now focus on the TS. The method of reduction of an infinite approximate Toeplitz 



matrix into a small finite one has been described elsewhere ||18|| . We thus only briefly sketch 
the solution. We have to substitute the exponential correction to the 3x3 matrix; the new 
matrix equation is : 
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(6) 



+ f h - 1 

where 6L and 9+ are complex exponents characterizing the (exponential) decay of the t n 's 
on both sides of the central energy band, n = 0), 9ft(0+) > 0). These exponents can 

be found by assuming that r n = r_i exp((n + 1)9-) for n < — 1 and r n = r\ exp(— (n — 1)9+) 
for n > 1, and substituting into Eq. @. We discuss three regimes : (a) fi 1 > 2^.fi (/c_i and 
k_ 2 are positive), (b) hQ < E < 2hQ (k-% positive and /c_2 imaginary), and (c) E < HQ 
(k_i and k_ 2 are imaginary). In the first regime the solution is : 



cos(^s(v±; 



(7) 



:)) = f (l + ^ + ^k\ 2 - ^(1-^ + ^ + ^(1+^)^) ^ 

cosh(K(0 ± )) = - £ cos(Q(9+))- 1 

while similar solutions can obtained for cases (b) and (c). 

It is clear from Eq. (|7|) that as e becomes smaller, the correspondence between FS and 
TS improves. Since we are interested in the limit e — >• 0, we regard the TS as an exact 
one. Fig. |2] compares the FS and TS in all three cases mentioned before (a-c) when e — 1, 
which is the worst case. There is a good agreement for the three central coefficients even 
in the extreme case of E jftVo- Significant relative deviations are obtained only for the 
levels whose population is exponentially small. Much better correspondence is achieved for 
smaller e, see Fig. 2d. Since in what follows we use only these three coefficients, this figure 
confirms the accuracy of the Toeplitz method for the following derivation. 
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B. Numerical Results 



We verify the above derivation by solving numerically the dynamical time dependent 
Schrodinger equation and comparing the results to the TS. In order to simplify the numerical 
calculation the delta function barrier is placed in the center of an infinite well. The algorithm 
we use is as follows. First, we find the eigenfunction and eigenvalues of the static problem, 
for each time. Once these are given, the numerical calculation becomes much simpler using 
adiabatic perturbation theory, leading to a set of first order differential equations. Due to 
symmetry, half of the matrix elements vanish. In the case of the delta function the odd 
eigenfunctions are time independent, thus the effective number of equations is reduced Q. 
The numerical solution is performed using the adaptive time step fifth-order Cash-Karp 
Runge-Kutta method [fI9|l . The initial wave function is a Gaussian wave packet located in 
the center of the left side of the barrier, moving toward the oscillating barrier. The wave 
packet collides with the barrier and splits to reflecting and transmitted parts. We measure 
the probability to be on the right side of the barrier after a collision. 

The same probability is calculated using FS in the following way. We average the trans- 
mitted probability current given by FS over the energy, using the weights obtained from the 
energy components of the initial wave packet taken above. The results of the analytic and 
numeric results are presented in Fig |^. The agreement between the two calculations is very 
good, confirming our derivation. 

III. TUNNELING TIME SCALES 

In order to investigate the tunneling time problem, one follows the strategy explained 
above, inspired by BL, and studies the frequency dependence of the transmission intensities 



The number of levels taken in the computation should include at least the first few energy bands 

(E) ± nhn. 
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in the limit of small e. In particular, one looks at the relative sideband asymmetry 

m) - Y^r 1 , (8) 

where I n = \t n \ 2 . BL show that for an opaque rectangular barrier, in the high frequency 
regime 

F{VL) = tanh( — — ) = tanh(fir BL ) (9) 

where tbl is defined as the tunneling time. The frequency = r^l is the transition 
frequency, below which the behavior is determined by the low frequency limit. 
In the low frequency limit, BL obtained the relation 

I ± = (eVr/2h 2 ) 2 (10) 

for an opaque finite barrier, whose height is V. Using this equation as a definition of r, they 
got the following result for a general barrier 
/ m 



T 



[k 2 - k 2 ) 2 n 2 d 2 + jfc*(l + K 2 d 2 ) sinh 2 nd + k 2 Kd(K 2 - k 2 ) sinh 2nd 1 1 2 



11) 



4fc 2 k 2 + fcg sinh nd 

where k = (2m/h 2 ) 1/2 VV - E, k = (2mV/h 2 ) 1 / 2 , k = (2mE/h 2 ) 1 / 2 and d is the barrier 
width. In the delta barrier limit this expression vanishes. The same result is obtained for 
the general rectangular barrier using the Larmor clock approach [11 



In what follows we use the first approach of BL, namely, looking at the asymmetry, and 
show that for a delta barrier in the deep tunneling regime E <g ^prV^ 2 , the frequency behavior 
depends on a new time scale. Using the above described Toeplitz method, we obtained for 
the asymmetry function F(Q) = —Qts where 

2h 3 . , 

T < = 7wf < 12) 

This holds only in the regimes (a-b), i.e., when E > Ml. In the other regime the behavior 



is also a function of ttrg, F(Q) = ^/Qts — To, where To is the tunneling probability in the 
absence of oscillations. Similar results for the model (0) were obtained by St0vneng and 
Hauge without using the Toeplitz approach 
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One thus sees that following the approach taken by BL, one finds a new and different 
time scale than the one obtained via the Larmor clock, which was also given in BL's original 
work using the low frequency limit. The interpretation of this result is not yet clear. We 
suggest that there are more than one time scale in the tunneling through a general barrier. 
BL claimed that since for the opaque barrier, in the high frequency regime one gets a 
clear definition of r = tbl (Eq. @)', this r is the tunneling time for the opaque barrier. 
Therefore, looking at the low frequency regime and knowing the tunneling time, one can 
extract a connection between the low frequency behavior and the tunneling time 

2h 2 i — 

r = %& (13) 



This relation was then used |13J to define the tunneling time in the general case. Now, if 
there is more that one time scale in the process, there is no reason to take the time scale 
which dominates the high frequency regime and connect it to the low frequency behavior 
which may be dominated by a different time scale. Looking at a delta barrier, in which the 
first time scale tbl vanishes and only the second time scale plays a role, one indeed sees that 
both the high and low frequency regimes are dominated by the same (second) time scale t§. 
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FIGURES 



FIG. 1. Oscillating delta function barrier. 

FIG. 2. Correspondence between FS (Full Solution) and TS (Toeplitz Solution), h = 2m = 1, 
Q = e = 1, V = 10. (a) £ = 2.5 > 2ft; (b) ft < E = 1.5 < 2ft; (c) £ = 0.5 < ft; (d) Same as (a), 
but for e = 0.5. 

FIG. 3. Tunneling probability as a function of ft - comparison of FS and dynamical computa- 
tion. The units taken are such that fi = 2m = 1, (E) = 5, Vo = 5, e = 0.9. The inset zooms on the 
maximum region 
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Fig 1 
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